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I. INTRODUCTION 

There has been considerable analysis [l| of the value of a high energy muon storage ring 
for the purpose of producing well-defined, high energy beams of neutrinos. This approach 
takes advantage of the relatively long muon lifetime to allow for collection at modest en- 
ergies and acceleration to very high energies, where the decay neutrinos form an intense, 
highly collimated beam. The proposals for such systems are quite expensive since the (gen- 
erally hadronically) produced muons must be captured, cooled and accelerated, all on a 
microsecond time scale. 

n 

Some time ago, Los Alamos proposed [2^ a pion post-accelerator (PILAC) to be appended 
to the then extant meson factory, LAMPF. The proposal took advantage of the high pion 
flux produced at modest energies, to collimate and inject into a high gradient linac. Phase 
rotation was introduced to reduce the energy spread of the beam considerably. The main 
practical difficulty with this relatively modest cost proposal was the short pion lifetime - 26 
nanoseconds. This required a very high accelerating gradient in order to achieve a useable 
flux at the exit of the acceleration section. 

The decay of a positively-charged pion is almost exclusively into anti-muons and muon- 
neutrinos (or muons and muon-antineutrinos for negatively-charged pions) with a purity at 
the part in ten thousand level. Thus a PILAC also provides an intense, collimated beam of 
neutrinos with an energy spectrum ranging from negligible to initially half the pion input 
energy and finally to half the output energy value. We show here that the neutrino spectrum, 
dNij/dE^, can adjusted to be almost flat in energy, if sufficiently high accelerating gradients 
for the pions can be achieved. The necessary values are below those already achieved 



a 



using superconducting RF cavities. [4] 

However, there is another advantage to such a neutrino source that does not dependent 
upon whether or not high gradients can be achieved. Currently, it is very difficult to de- 
termine the energy of the neutrino involved in a particular scattering event in the detector. 
Usually, a complicated series of Monte Carlo simulations of the detector and its response 
defines a statistical determination only. In the scheme we are discussing, however, it is 
possible to obtain the incident neutrino energy directly, event by event. 

We analyze here the theoretically simple case of a uniform accelerating gradient analyti- 



cally. A more realistic case with drift sections is left for the accelerator physicists who would 
actually design and build such a system. 

II. CALCULATIONS 

Although we ultimately want to obtain the neutrino energy spectrum as a function of 
distance from the injection point of the accelerator, to follow the decay rate, it turns out 
to be more convenient to first find the pion energy and decay rate as functions of time. 
Throughout, we will use units where h = c = 1. 

A. Time Dependence 

We begin with pions injected at energy Ej and 7/ = Ej/itIt^ where 771,^ is the pion mass. 
For simplicity, we assume a spatially constant, uniform accelerating gradient, G, even though 
the actual acceleration structure will consist of short, high acceleration regions interspersed 
with drift regions. Refinement to include this detail should only change our smooth curves 
to segments of alternating flat and higher slope regions but with essentially the same average 
behavior. For convenience, we scale G by the pion mass so that g = G/itLj^ has the dimensions 
of an inverse length, 

grriTr = dE^^/dx . (1) 

The energy of a pion as a function of distance along the accelerator is therefore 

E^{x) = Ei + gm^x . (2) 

However, to follow the decay process, we will need the time rate of change of pion energy, 
viz., 

dE^/dt = -7^-77 = 9^-Kl3{t) , (3) 

where (3{t) is the pion velocity in units of c. Since 7(t) = ET^/niT^, we can rewrite Eq. (1) as 

d^/dt = g(3 (4) 

and since 7 = 1/Vl ~ P'^, 

d^ ^ (3 d(3 

dt [1 - p^] 3/2 dt ' ^ ^ 



Eq. (4) thus becomes 

This is easily integrated to obtain 



1 dp 



[1 - /32] 3/2 dt 



\m=9t + 
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or equivalently, 



This allows us to identify 



[l-/32(t)]l/2 ^ [1-/52(0)]V2 ' 

/3(t)7(t) = ^t + /9(0)7(0) = (7t + «o • 
gt + ao 



(3{t) = ■ 

^l + igt + a^y 

7(t) = ^l + ^gt + aoY 



(6) 

(7) 
(8) 

(9) 
(10) 



Bear in mind that ao is implicitly a function of the injected pion energy. Eq. ([8]) is just 
the well-known result that a constant accelerating energy gradient causes the momentum to 
increase linearly with time, i.e., p{t) = iriT^gt + p/, where pi = mjrCto is the pion momentum 
at injection into the accelerator. 



B. Decay 

With 7(t) in hand, we can solve for the pion decay function in the rest frame of the 
accelerator, viz. 



dN^ 



r 



■iv.(t) 



-N^{t) 



(11) 



(^t lit) ^i + ^gt + aoY 

where F is the decay rate of the pion in its own rest frame. Integration of this equation 



yields 



N.{t) = iV^(O) 



(gt + ao) + V 1 + (gt + ao) 



-r/5 



(12) 



ao + y 1 + "o 

Since the number of neutrinos is just equal to the number of pions that have decayed, 
N^(t) = iV^(O) - iV^(t), it follows immediately that 

r/g\ 



N,{t) = iV^(O) 1 



ao + V^ + '^o 



{gt + ao) + Vl + {gt + «o)^ 



(13) 



C. Distance Dependence 



We next convert this time dependence into spatial dependence, as a function of distance 
from the input of the hnear accelerator. We can integrate (3 = dx/dt given by Eq. (Q to 
obtain 



x{t) = / dt'/3{t') 



1 + {gt + ao) 



1 + a^ 



/o g 

where we chose x{0) = 0. Eq. flT4|) can easily be inverted to obtain 



tix) 



gc 



(14) 
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and this may be substituted into Eq. flT3l) to obtain Ny{x). Our result for the neutrino flux 

at a distance x is, therefore, 

r/51 



N,{x) = iV,(0) 



gx+Jl + al + [{gx + Jl + a2 )2 _ 1 ] 1/2 



Noting that m-,^Jl + Oq = Ei^we convert this back to the original variables. 
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K{x) = NM 



NM 



Pi + Ei 



r/g 



1 



^gm^x + Ej + ^[gmT,x + E/]^ 

- J-, , s mT/G-\ 

Ej+pi\ 



mt 



(17) 



E^+Px) 

where E^ and p^ are the energy and momentum of the pion after travelling a distance x 
through the accelerator. 

Fig. [1] shows the behavior of N^{x) for two different accelerating gradients and for two 
different initial pion kinetic energies (in the laboratory). As expected, the x-dependence 
becomes more and more linear as the injection energy increases. Also, the higher the gra- 
dient, the flatter the curves become. The reason for this is that a bigger gradient means 
higher energy pions, which have longer lifetimes (in the laboratory), whence fewer decays 
into neutrinos during acceleration. 

The accelerator designer may therefore want to optimize between the number of neutrinos 
showing up in the detector (normalized by the number of pions captured in the accelerator) 
and their ultimate energies. One might also want to consider the possibility of making a 
pion (ie-celerator, if it turns out that lower energy neutrinos have desirable properties for 
some experiments. 



D. Neutrino Energy Spectrum and Angular Divergence 

The neutrino energies will range over a wide spectrum, but for any specific time from in- 
jection of the pion bunch and angle of divergence from the beam axis of a neutrino scattering 
event in the detector, only one specific neutrino energy is possible. 

The momentum and energy, Eq, of a neutrino emitted in the pion rest frame is 

2 2 

Eo = -^ ^ = 29.8 MeV . (18) 

2m^ 

For a decay where the neutrino is emitted at a polar angle Or (in the rest frame) with respect 
to the pion beam axis, the momentum, polar angle 6l, and energy E^, in the lab frame are 
given by 

p = EoW + cosOr), sin^^J , (19) 

tanOL = smeR/l-fiP + cosOR}] , (20) 



E, = Eo^sm^eR + ^^i/3 + coseR)^ 

= -f{l+i3coseR)Eo, (21) 

where we show only the longitudinal and transverse momentum components, ignoring the 
azimuthal angle since the distribution in that angle is uniform. 

From Eq. fl?Il) . neutrinos emitted exactly backward from the beam direction will have lab 
energies and momenta of 7(1 — /5) Eq, i.e., approximately zero for the most relativistic pions. 
Those emitted precisely forward will have lab energies and momenta of 7(1 + /5) Eq. For 
neutrinos emitted transverse (in the pion rest frame) to the beam direction, the lab energy 
and momentum will be '^Eq. 

From Eq. (BU]) . the singularity at cos6r = —(3 corresponds to a laboratory angle of 90°. 
Neutrinos at this angle come from a very backward decay and have a low energy of -Eq/t- 
Those that come from even more backward decays have yet lower energies. The maximum 
laboratory angle is 180°. That is, 9^ as a function of Or rises monotonically from to 180°, 
and does so more steeply in the backward Or direction. 

Neutrinos are emitted between the two extreme energies along a curve in Ey-O^ space. To 
make such a plot, one needs to find the value of cos Or corresponding to a given laboratory 



angle, 9l. From the inverse Lorentz transformation, from the laboratory to the pion rest 
frame, and with some algebra, 

cos Ol — P 



cos 9ji = cos 



tan ' 



7(cos6'l - /3) 



1 — /3 cos 6l 



(22) 



Having found the value of cos Or for a given (or derived) 6l, the neutrino energy E^, is then 
given by Eq. (pTj) . or more directly, E^, = (1 — /?cos6'i)^^-Eo/7- Figure [2] shows curves of the 
neutrino energy as a function of 6l for three different (total) energies of the pion at the time 
of its decay. Note that, for the higher pion energies, the neutrino energy Ei, is more forward 
peaked and, of course, is larger there. Taking into account that neutrino cross sections are 
proportional to Ej^ at low E,^, this suggests that, for economy, the neutrino detector could 
well be long and narrow, subtending a small solid angle along the beam direction. 

E. Determining E^ From Event Timing 

Fig. 3 depicts an example of what we have in mind, a pion accelerator which is 10 m 
long. [5] This is followed by an air- filled drift space (i.e., a decay volume) of 14 m, a 10 
m long beam dump, such as earth or iron, and the neutrino detector, presumed, for cost- 
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effectiveness, to be a long cylinder aligned with the beam axis. J6|] The drift space could also 
provide space for experiments using the pions which emerge from the accelerator. It also 
might contain bending magnets to deflect charged particles from the beam axis. [3] 

Neutrinos all travel at velocity c, which is discernably larger than the velocity of the 
pions (up to ^ 1 GeV pion energy). Thus the time of scattering and the position of the 
event in the detector are uniquely correlated (up to resolution limits) with the energy of the 
neutrino. However, low-energy neutrino cross sections vary as E"^, so most of the neutrino 
events will come at forward angles less than ?a 30°. 

We discuss first the case for those higher-energy neutrinos emitted at 0° along the beam 
direction. Fig. [3] shows how the change in pion velocity along the accelerator, followed by 
speed-of-light propagation of the neutrino from the decay point, determines the neutrino 
energy. Note that this figure is for an unrealistic injection energy of T^ = 100 keV, simply to 
show the principle involved. [Sj] The idea is to measure the difference between the time when 
the neutrino scattering event occurs and that for a light-like signal that begins at the time 



the pions are injected into the pion accelerator. The neutrino energy E^, in this on-axis case 
is given by 



E, 



u, on— axis 



7(tdec) [1 + /?(tdec)] Eo 



(23) 



where tdec = t{xdcc) is given by Eq. (TTS!) and Eq by Eq. (TTSj) . The value of Xdec needed for 
tdcc is obtained by solving the more general expression for the event time, Eq. flMl) . given 
below. 

Note also that some pions will decay after leaving the accelerator. In such cases, as in 
the dot-dashed curve in Fig. 3, the timing for the event will be later than the case shown 
when Xdecay = 10 m, but E^, will be the same, i.e., 93.0 MeV in this example. 

For more realistic injection energies, e.g., T^ = 350 MeV, the curves in the pion accelerator 
section are much, much flatter than those shown in Fig. [3l The time differences are quite 
small but technically feasible. For an injection energy of T^ = 350 MeV, Table I shows, for 
two choices of accelerating gradient G, the time differences, in nanoseconds, and the higher- 
energy neutrinos' associated energies. Note that, to distinguish the two highest energies 
shown for G = 50 MeV/m, one must have timing resolution of about 50 picoseconds. 

A similar but more complicated relation occurs for radially off-axis points, where the 
decay angle must be included in the correlation. (See Fig. [2l) The procedure for determining 
El, from the time the event occurs (relative to the time of injection of the pions into the 
accelerator) is as follows. 



^decay i^^j 


G = 10 


MeV/m 


G = 50 


MeV/m 




At (nsec) 


E^ (MeV) 


At (nsec) 


E^ (MeV) 


0.0 


0.000 


205 


0.000 


205 


2.0 


0.275 


214 


0.236 


249 


4.0 


0.528 


223 


0.401 


292 


6.0 


0.762 


231 


0.522 


335 


8.0 


0.978 


240 


0.616 


379 


10.0 


1.179 


249 


0.691 


422 



TABLE I: The calculations shown here are for T^r = 350 MeV at injection into PILAC. 



If the position of the event in the detector is given by (xcvt, Pevt), where x^vt is measured 
from the injection point and pevt is the transverse distance from the beam axis, the time of 
the event is 

^evti-^dec; -^evt) Pevt) ''dccl.'^dccj T i^u[,-^evt -^deo Pevt) ; \ / 

where tdec(a;) is given by Eq. IHM and t^{x,p) = y/x'^ + p'^/c. Equation (^^ can be solved 
numerically for Xdec- The implicit dependence in tdec on the gradient, G, and the kinetic 
energy T^ of the injected pions should be kept in mind. 

Given the extracted value for x^ec, we can infer that the neutrino was emitted at a 
laboratory angle of 

9l = tan~^[pevt/(a;evt - a;dec)] ■ (25) 

The pion energy at the time it decays is also determined by Xdec- From Eqs. ( TTOl) and (TTSll . 



7dec = \/l + \gt{Xdec)+C^ , (26) 



/?dec = V^dec - 1 /7dec , (27) 

Ejr = m^ 7dec • (28) 

From the values of 9l, /^dec, and 7dec the cosine of the decay angle in the pion rest frame 
is then found by substituting these values into Eq. ( l22l) : 

COS^L-Z^dec .r>n\ 

1 - /idee COS Ol 

The neutrino energy E^, is then determined by substituting this value of cos9r into Eq. fl2Tl) . 

Table II shows the correlations between measured scattering times At and off-axis event 
positions pevt and the extracted values for the decay point Xdec, laboratory neutrino angle 6l, 
and E,^, for the same two accelerating gradients and injection energy used in Table I. Because 
of the small solid angle subtended by the detector, in this example, the time differences 
at off-axis points are not much different from those shown in Table I. Thus transverse 
spatial resolution becomes a significant contributor to determination of the neutrino energy, 
implicitly demanding a fine-grained detector. 

Fig. 4 illustrates the correlation for an off-axis detection of a neutrino between the time 
of the event and the pion decay point and angle of decay (as seen in the lab frame), whence 
the energy of the reacting neutrino. 
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At (nsec) 


Pevt 
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3^ dec 


Ol 


Eu 


0.250 
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1.81 
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213 


2.15 
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252 




0.6 


1.72 


15.7 


213 
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15.8 


249 
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1.58 


23.4 
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1.82 
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245 
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3.77 


8.3 


222 


5.58 


8.7 


327 




0.6 


3.67 


16.5 


221 


5.35 


17.3 


322 
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3.50 


24.7 


221 


4.99 


25.7 


314 


0.750 


0.3 


5.89 


8.9 


231 


11.9 


10.7 


463 




0.6 


5.77 


17.5 


231 


11.3 


20.9 


450 
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5.58 


26.1 


230 


10.5 


30.6 


433 



TABLE II: Extractions of the decay point x^ec (in ni), neutrino laboratory angle 9l (in mrad), 
and El, (in MeV) for a scattering event taking place at various scattering times At in the detector 
at three values of off-axis radial position Povt- This case is for an event occurring in the detector 
at Xevt = 40.0 m from the accelerator injection point. We have assumed here that Tt^ = 350 MeV 
at injection into PILAC, as in Table I. 

III. DISCUSSION 



Note that, as shown by Eq. fITTj) . the important scale for the acceleration gradient is set 
by the lifetime of the pion in units of pion mass, namely. 



G = m^T^/c = 17.9 MeV/m . 



(30) 



Gradients exceeding this value are now routinely available, which makes possible very tightly 
focussed neutrino beams and unusual neutrino energy spectral distributions. 

The higher the gradient can be made, the more pions survive to the exit energy and the 
neutrino spectrum becomes more and more weighted to higher energies. The higher energy 
component of the neutrino beam also becomes more focused, projecting a smaller image on 
the target. 

The pion particle flux expected from the design in Ref. [2.] was lO^s^^ at the target end 



of the accelerator. With the gradient given above, the accelerator would be on the order of 
5 lifetimes in length. Thus, while the neutrino flux at the highest energy available would 
be on the order of 10^s~^, the total flux over the entire spectrum would be on the order of 
lO^^s"^ (since nearly every injected pion will decay into a neutrino). This certainly presents 
a challenge to experimenters to design detectors with a useful event rate, as the flux is not 
as large as produced by simple beam target-decay volume sources. 

An additional order of magnitude in flux should be available if the PILAC system were 
to be mounted at the Spallation Neutron Source at Oak Ridge National Laboratory. 

We thank J. Jenkins, W. Louis, G. Mills, and C Morris for valuable conversations. This 
work was carried out under the auspices of the National Nuclear Security Administration 
of the U.S. Department of Energy at Los Alamos National Laboratory under Contract No. 
DE-AC52-06NA25396 
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FIG. 1: Dependence of N^{x), normalized by the number of input pions, Nj^{0), as a function of 
distance along the pion accelerator, x, for two pion kinetic energies spanning the production peak 
at injection into the accelerator. The solid curves are for a gradient of 10 MeV/m, while the dashed 
curves are for 50 MeV/m. 
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FIG. 2: The neutrino energy, E^, (in MeV), as a function of laboratory polar angle 0l from decay 
in the acceleration section at a given time from bunch injection (or position in the accelerator) for 
three pion kinetic energies at the time of decay. Each pion energy corresponds to a specific location 
for the decay in the accelerator section. 
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FIG. 3: Relation between decay time after pion bunch injection and time of scattering event in 
detector for on-axis u^s of the highest energy. This is for an unrealistic injection energy of T^^ = 
100 keV, so that the curvatures before decay are highly exaggerated. 
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FIG. 4: Relation between the reacting neutrino's energy, Ey, and the time difference between the 
time of bunch injection and the scattering event time in the detector for a particular event position, 
here at 40 meters from the start of the pion accelerator and at an off-axis distance of 1 meter. It 
shows the differences in Ey and time difference At for three different pion decay positions. The 
kinetic energy of pions at injection has been taken to be 350 MeV and the accelerating gradient is 
10 MeV/m. Thus, in this example, if the absolute time of this scattering event can be determined 
relative to the injection time with an accuracy of better than ~ 500 ps, then the energy of the 
incident neutrino that initiated the event will be determined to an accuracy of better than ~ 15 
MeV. 



